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STROP SIS 


Fading dispersive channels, like troposcatter/ 
iono scatter channels, mobile to mobile platform communication 
channel via satellite and satellite to satellite relay channels 
are in general characterized by a linear time-varying filter. 
Corresponding to the transmission of a sinusoid over such a 
channel, the received signal contains a specular component and 
a random component which is characterized by a nonstationary 
narrow band Gaussian process. Over short time intervals the 
random process may be approximated by a stationary process which 
has a Rayleigh distributed envelope. Hence the composite signal 
will have a Rician distribution. To combat fading the same 
signal is transmitted over several paths and various diversity 
combining techniques are used at the receiver. In general, 
both the specular and diffused signal components may be 
different over various diversity paths. For example, in angle 
diversity troposcatter system the diffused component in the 
signal corresponding to the antenna beam at grazing horizon is 
maximum and is minimum for the beam with the highest elevation. 

In this thesis for digital communication over such 
channels, the expressions for the probability of error are 
obtained asstiming post detection diversity combining. A general 
type quadratic form detector is considered which includes 
incoherent, differentially coherent and incoherent 
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detectors as special cases. The probability of error expressions 
are evaluated by applying the method of residues. The method 
consists of finding the two sided Laplace transform, F(s), of 
the probability density function of the sufficient statistic. 

The residues of F(s) in the right half complex plane are 
evaluated to yeild a function p“ (w) that vanishes for w;>o 
and whose integral gives the probabilities of error. In those 
cases where F(s) can be easily split into two parts F^(s) and 
F^ (s), (as for example when F(s) is a rational function) which 
are analytic in the right half and left half s-plane respectively 
then even the calculation of residues is dispensed with and the 
probability of error is simply obtained by evaluating F^(s) 
at s=o. Due to the more direct nature of method, the resulting 
expressions are considerably simpler than the ones obtained in 
the earlier literature for the situation where the diffused 
components are equal over various diversity paths. The expression 
derived consists of a single convergent series of which only 
a few terms need to be evaluated in most cases of practical 
interest. Further Hight upper and lower bounds are found on 
the series sum which consist of a sum of a few polynomial, terms 
of finite degree for the complete range of signal to noise ratio. 
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The result of the generalized quadratic form 

receiver are specialized to slow fading Rician and Rayleigh 

channels for FSK, DISK and coherent PSK modulations. High 

SHR assymptotes are found in each of these cases. The 
of the 

expressions/ prob. of error obtained are either same as 
derived earlier in the literature or are considerably simpler 
in form and easier to compute. The final expressions do not 
involve any special functions like generalized r function and 
thus require no special computer sub-routines [39]. 

In this thesis also is analyzed a more general 
situation in which the total number of diversity channels could 
be divided into F groups such that the channels within a group 
have equal diffused components. This model for example may 
correspond to dual- space dual angle diversity troposcatter 
system. The expression is simplified further when no two 
channels have equal diffused components as may be the case in 
triple angle diversity troposcatter system. 



CHAPTER 1 


INTRODUCTION 

1.1 A MODEL FOR FADING- DISPERSIVE CHANNELS s 

Fading dispersive channels like tropo scatter ./ 
ionoscatter channel, mobile to mobile plateform communication 
channel via satellite and satellite to satellite relay channels 
are in general characterized by a linear random time-varying 
filter. Corresponding to the transmission of a sinusoid over 
such a channel, the received signal contains a specular 
component and a random component which is characteriz ed by a 
fconstationary narrow-band Gaussian process. Over short time 
intervals the random process may be approximated by a stationary 
process [l-7] which has a Rayleigh distributed envelope. Hence 
the composite signal envelope will have a Rician distribution. 

The deterministic and random components may arise 
due to different sources in various physical channels. For 
example the random component in tropo scatter channel is 

considered to occur as a result of scattering from 
a large number of scatterers in the common volume between the 
transmitter and receiver antenna beams. The specular component 
in these channels may be due to ducting phenomenon which may 
normally occur in tropical or subtropical climates [8-9]. In 
mobile communication via satellite, as for example in 
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communication from an aircraft to another aircraft or to the 
ground station, the signal has two paths. The first path is a 
direct LOS path from aeroplane to satellite which has fixed 
mode propagation and results in the deterministic component 
of the signal. The other path consists of signal reflected 
from the ground and then reaching the satellite. The nature of 
this signal in most cases is diffuse especially when the 
terrain below the aircraft is rough. As shown in [10-12] this 
component has a Rayleigh distribution. Hence all these channels 
are modeled as shown in the figure 1. The block entitled 
deterministic filter corresponds to fixed mode propagation 
path, while random filter corresponds to fading mode propagation 
In the figure s(t) is the complex envelope of input signal and 
r^ (t) and r c ^ (t) are complex envelopes of fixed mode or 
specular component and diffused component of the received signal 
over k th diversity branch. At the receiver complex gaussian 
noise is added to give the complex envelope r^.(t) which 
constitutes input to the quadratic form receiver. 

To combat fading, same signal is transmitted over 
several paths and the various diversity combining techniques 
are used at the receiver. In general both the specular and 
diffused signal components may be different over the various 
diversity paths. For example in angle diversity troposcatter 
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system the diffused component in the signal corresponding to 
the antenna beam at grazing horizon is maximum and is minimum 
for the beam with the highest elevation. 

1.2 general Quadratic Form Receiver : 

Figure 2 shows the block diagram of the quadratic 
form receiver [13] for which the probabilities of error are 
evaluated when operating over Rician fading channel. The 
outputs U k (t) and V k (t) are complex valued Gaussian processes 
obtained as a result of linear operations upon r{t). The 
processes U k (t) and J k (t) could be matched filter outputs, or 

U k (t) could be a filter output and V k (t) a reference waveform. 
The outputs 2 k (t), Y k (t) are passed into a device which 
computes the quadratic form 

w k (t) = f |n k (t) | 2 + g |Y k (t)| 2 + eU* it)V k (t) 

+e U k (t) I k (t) (1) 

where f and g are real and e may be complex. 

The quadratic form (l) is summed, together with 
corresponding quadratic form contributions from other diversity 
channel to yeild 

w(t) =5Tw k (t) (2) 

w(t) is sampled at time instants t=mT(m being an 
integer and T is symbol duration) and is compared with i-ero. 
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FIG. 2. QUADRATIC FORM RECEIVER 
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The transmitted symbol is decided to be ’"l” if w(mT)^>0 and 
" 0” otherwise. 

When the channel exhibits frequency selective fading 
for high data rates there will in general be intersymbol 
interference from adjacent symbols. One of the possible ways 
to take into account the intersymbol interference is to 
calculate the probabability of error for a fixed pattern of 
interferring symbols and then average over all possible patterns. 

let P-^ q ( ^ ) denote the probability of deciding ! 

" 0 M when actually *' 1" has been transmitted conditional on 
the digital sequence which contains ” 1” in the 

particular decision interval. 

Then the average probability of error is given by 

Pe = | P. +| 21 p 01( S^.^ 0 ^ )p ( S-^.^ 0 ^ ) (3) 

where is a digital sequence having H 0 1 ' in the 

decision interval 

Pq^ is the probability of error when M 0 " has 

been transmitted. 

p(S, (l) ) is the probability that the transmitted 

- « ( 1 ) 

sequence is S, . 

By finding and V^. for each of possible sequence 

or one can calculate P-^q (S^^^or ^Ql^k^^ an ^ 
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then from equation (3), P@ could he obtained. 

1.3 Brief Review of Prob. of Error Calculations : 

Probability of error analysis for slowly fading 
multipath and diversity channels having no specular component 
have been carried out by Pierce [14] , Turin [15]. Turin [16] 
considers an optimum diversity reception through dependent 
Rayleigh fading paths with no explicit channel measurements. 
Pierce and Stein [ IT ] consider the analysis for optimum 
predetection combining when the channel undergoes dependent 
Rayleigh fading that is presumed to be measured perfectly. For 
independent fading, Bello and Relin [18] derive expressions 
when the measurements are noisy. Further, Proakis, Drouilhet 
and Price [19] give experimental results pertaining to the 
performance of coherent detection system using decision 
directed channel measurement. Lindsey computes the error rate 
for coherent and incoherent systems operating through slow 
fading Rician multichannel. 

Price [20] considers the fast fading scatter channel. 
Bello and Relin [21] have examined matched filter receivers 
subjected to fast fading and have shown that an irreducible 
error probability exists over such fast fading channel. Pierce 
[22] has computed error probabilities for certain spread 
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channels and Kennedy [l] has developed certain performance 
hounds over some scatter channels. More recently Monsen [23], ' 

[24] has developed adaptive equalizer for such fading channels 
and has computed error rate performance and has shown an 
implicit diversity improvement due to dispersion. 

Most of this literature is however confined to channels 
which don't have any fixed mode component. Lindsey [25-28] and 
Bello [13] have derived expressions for probability of error over 
Rician dispersive channel when all the diffused components are 
equal . 

However there are physical situations, for example, 
angle diversity troposcatter systems [29,30], where the 
diffused components received over various diversity paths are 
•unequal. 

1.4 Outline of the T he sis s 

In this thesis expressions for the probability of 
error are derived for digital communication over such channels, 
assuming post detection diversity combining. A general type 

quadratic form detector [13] is considered which includes 
incoherent, differentially coherent and incoherent detectors 
as special cases. 

In this thesis a different approach is employed to 
calculate the probability of error expression over dispersive 



9 


Rician fading multichannel , which results in expressions which 
are considerably simpler than those derived by Lindsey [27] 
and Bello [13]. Further the method admits the generalization 
to the situation when both the specular and diffused components 
are unequal over different diversity paths. The method consists 
of finding the two-sided Laplace Transform F(s) of the 
probability density function of the sufficient statistic. The 
probability of error Pe is then calculated by finding the 
residues of F(s) in the right half complex plane getting a 
function p~(w) which is zero for w;>0 and then integrating 
p"(w). In those cases where F(s) can be easily split into two 
parts F-^s) an & F^Cs) a ® ^ or example when F(s) is a rational 
function, which are analytic in the right half and left half 
of the complex s-plane respectively, then even the calculation 
of residues is dispensed with and the probability of error is 
simply obtained by evaluating F^s) at s=0. This method differs 
from the earlier approaches in recognizing that to calculate 
Pe it suffices to find p~(w) from F(s) rather than p(w). 

In chapter 2 the probability of error expression is 
found for the general quadratic receiver. From this the 
probability of error for selective fading channel could be 
found by calculating Pe for each possible sequence of interferr 
ing symbols and averaging the result. For the general quadratic 
receiver an upper and a lower bound is found which consists of 
L number of polynomical terms, each having a degree less than L 
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In chapter 3 various results for the slow fading 
Rician and Rayleigh channels are derived as special cases of 
the more general expression derived in chapter 2. Both coherent 
and incoherent detection schemes are considered. The specific 
results agree with the ones derived earlier in the literature. 

In chapter 4 a more general diversity situation is 
considered where the diffused components are different over 
various paths. Prohahility of error expressions are derived 
for a slow fading Rician multichannel configuration with 
incoherent and partially coherent reception. The situation 
considered is the one in which the total number 1 of diversity 
channels could be divided into N groups such that the channels 
within a group have equal diffused component. This model for 
example may correspond to dual space-dual angle diversity 
troposcatter system. The expression is simplified further when 
no two channels have equal diffused component as may be the 
case in triple angle diversity troposcatter system. The result 
of chapter 3 is shown- to be special case of this general 
result. The probability of error expression is also found for 
selective fading Rayleigh channels with unequal components and 
employing general quadratic form receiver. 

Chapter 5 concludes the thesis with some suggestions 


for future work 
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CHAPTER 2 


PROBABILITY OP ERROR WITH QUADRATIC P0BM RECEIVER 

2*1 In this- -chapter. expressions' for the probability of 
error are derived for the quadratic' -form receiver considered, in 
section 1*2 A general expression is 'derived which. consists of 
.a convergent series which, converges faster, than an exponential 
series. Por most of the cases of practical. interest only a few 
terms of the series will.be required to "yeilcL the pro.ba bi lity 
of error. The general Expression is bounded by an upper and a 
lower bound each of which contains summation of finite number 
of terms. 

2,2 General Expression for Probability of Error » 

As mentioned in section 1.2 to calculate probability of 
error for frequency selective fading" channel from equation (3)» 
one has to calculate the probability of error and 

Pqi( ) for each possible interferring sequence S 

However, in the following we will derive P^q (-0^^') or Pq^ ‘ 

for a particular sequence so that the statistics of the 
resulting signal in the detection period is known. Averaging 
this error over the possible sequences will yeild the. required 
probability of error. Por simplicity of notation P 10 
or P 0 i(^°)) which will be equal under equal signal energy 
assumption will simply be denoted by Pe. Assuming the 
equiprobable binary signals, the threshold of the receiver 
will be at zero. 
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How denoting w(mT) by w , the probability of error 

is given by 

Pe = Pr [ w<0] 

0 

= I" p(w) dw (4-) 

Using Turin's resylts [31] » Bello [13] has derived 
the characteristic function of w^.. Here, we will use the two-, 
sided laplace transform of p(w k ) instead which is given after 
a slight change in the form by 

p k a k 

exp c i +®k~ ] exp [ ' r ^; ] 

-«*rtp(* k )]- K k — ^ (5) 


TVfap- 


V d k } 


where K fc = a fe d fc exp (~P k /a k ) exp (- 

a k and d k are evaluated as solution of the following 
equation 

|l- sH^G j (s+a k ) (-s+d k ) 


( 6 ) 


and p k and q k are calculated in terms of first and second 
order moments of U k and V^as follows. 


Letting G- = 

f e" 

- e S'. 

ii 

M 

W 

n u m uv" 




m uv m v ; 


and S k = 

where f and g are real and e may be complex. 


Hi 


T?i 
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and where 


n U ” 1 2k “^k 


m ? " -tf 


n ui = ( y k *k> 


^k ~ -K » **k ~ -k 


one computes the matrices 

2d 


G 


ak 


G 


A 


hk 


( 


k 


d k +a k 


2S.n 


. d k + a k 


)• ( H, 1 - a, G ) 


) ( ^ ) 


( 7 ) 


then the p k and q k are given by 

1 


p k = 


2 a k S k ®bk S k 


q k = | d k S k G ak S k 

and from equation (2) one obtains 
I 

F (s)= TT P k (s) 

k=l 

Now let F(s) be split such that 
P(s) = F-^s) + F 2 (s) 

f!(s)=/ [ P + M ] 

t p 3 


( 8 ) 


( 9 ) 


where 



U 


and 

P + (w) = p(w) w ^ 0 

= 0 w <. 0 

P ~(w) = p(w) w 0 

= 0 w > 0 

The function p(w) being the probability density 
function which is a convolution of several continuous bounded 
functions, is bounded over the entire range of w. Hence 
F-^(s) and T^ 13 ) cannot have any singularity in the right and 
left half of s-plane respectively. 

As p + (w) and p“(w) together constitute p(w), we have 
p"(w)= -sum of residues of P(s) at the singularities in the 

right half of the s-plane (10) 

p + (w)= sum of residues of P(s) at the singularities in the 

left half of the s-plane. 

In the situation where all the diversity channels 
have equal diffused components but different specular components, 
each of the diversity branches will have the same H-^matrices 
because nig, m-y and nigy don’t depend upon the means p-^. and ^ ^ 
which will depend upon the specular component. 

Accordingly one obtains 

a^ a , d k =d for k=l,2...,L 

Hence from equation (9) after substituting (5)»one obtains 

Is 1 C rsV- ^ 

P(s) = K 


( 11 ) 


(s+a) L (-s-H) L 
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where E = 


and p = 


7“ K, = a JJ d P exp (-p/a) exp(-q/d) 

4rh ♦ 

k=l 

1 L 

fc-, P k ' q = SI q k 

x k=l 


Now the singularities of E(s) lie at s=-a and s=d 
where a and d are positive real numbers. 

Hence from equation (10 ) one gets 

p“(w) =- residue of F(s) at s=d 
or p“(w) = coefficient of TLs'-fc l ) in 'the Laurent 

series expansion of F(s) e sw (12) 


Expanding exponential terms in equation (11) one obtains 

a> 


P(s) e sw =K 21 — i — - -JL. 

A v to / ^ -to- -! — n TT.J-t -5 t 


1=0 (OtaT I+i i! [ Jtt-a) 1 + (d- s ) 1+1 + +, "‘ ] 


2 

X-- C l-t(d-s)(-w)+(d-s) 2 i^j— + ,...] e 


wd 


(13) 


Now letting M=L+i and repeatedly using the expansion 


s 45-= (4+arl [ Hni r 1 

= (a+a)- 1 [ i+ -!=!- +.... ] 

in (13)? one obtains. 


(14) 


E(s) 6®"= E [ 


CQ 


P 


M, . f 1+P M,l ( ' d+a ^ +P M,2^ d+a ^ +,# *^ 


X- i 


Xiq (d+a) i! 

1 __1_ 

(d-s) 1 + (d-s) L+1 " 1 


r4r +•••.. ] 


X' [ l+(d-s) (--w)+(d-s) 2 ( -w .) 2 + 

2 ! 


• • * • 


] e wd ] 


(15) 



(M+n-l) ! 

where P.. = — 7 ^ v , - — 

M,m (M-l ) !n 1 


for M>1, m>0 and P^ m =0 for m < 0. 


Multiplying out the first two series and arranging the result 
in powers of (d-s), one obtains 


F(s) e sw = E X 


co 


i =0 (d+a) 


M 


[ H (a-s)I+ xiri) [ 


'M.l-1 


i! 


1 1! 


x M.l 
(d+a) 

1 

(d-s ) 2 L (d+a ) L ~ 2 


r M . 1+1 
(d+a) 


?=0 


. 2 
.•a_ 

L+l 2! 


r P M,L-2 , P M.1-1 

L TO* 


(d+a) 


1-1 IT 


“tT 


~M,1 


(d+a)' 


(d+a) 


1-1 


-§r+ ...] 


(a- 


+ 


(d-s)' 

1 

(d-s) : 


p p 

r M.l . . M.2 

ci , 

P M t 3 

^ d+a} (d+a) 2 

Jr + 

(d+a) 

, P M,1 Q_ . 

p 

r M.2 

2 

„ q 

[,0 + Td+a) 1! + 

(d-m) 2 

2! 

P , p M t i 

2 

cl. . 

P M«2 

P M,0 1! + (d+a) 

2! + 

(d+a) P 

2 ( — W ) 2 



r) + (d-s) 

Hf" * m • * 

] e 


3 2! 


.. 3 


3 


7! T * ♦ 


wd 


• 3 +• • • 

(16) 


The precise values of in ( 16 ) are of no interest 
in what follows. 

The coefficient of rg— 7 ) appearing in the equation 
( 16 ) which is also equal to p*"(w) is given by 
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p-(w)=K fb -At- [ [ h + 1%.-] 

fe) (d+a) M i! (d+a) L “ 1 (d+a) 1 1 * (d+a) 1+1 2 * 

. (-w) r P H:i-2 . P M.L~1 a P M.I a 2 . n 

+ — * — Y~o + — * — T~1 n + t" oT +• . . J 

11 (d4a) L “ 2 (d-fa) 1 " 1 1 - (d+a) L 2 ' 

+ LSL^r 5Lh=2 + ■ P M,L-2 a P M , L--1 4 

2! [ (d+a) 1 *^ (d+a) 1 - 2 57 (daa)^ 1 ** *** ] 


+ 


+ _lzw4l 1 r P + . f M y i 

* * * + (1-1)1 L X M,0 + d+a 


SL. + JLI 

- 11 ' ' 
- 1 - (d+a)' 


4 i + ... .3 


+ 


L-yJ: 

Ll 


r p 2— + —JLti. 

L X M,0 1! + (d+a) 


4 + ^ z 

(d+a) 2 


-^T+. . . . ] 


+. 


] e 


wd 


(17) 


Using the identity 


f '4r eWd d " 


a 


1 _ 

X4-1 


(18) 


One gets after substituting for p" - (w) in equation (4) 


Pe=K 


? P r 1 r P M.L-1 . P M.L a ^ P M.L+1 q 2 ^ 

iS5( a+ S) M i! [ d [ (d-m) 1 - 1 (d^) L ^ (d +a ) I+1 21 ] 


d' 


r P M,L-2 . P M , 1-1 a . P M,L q 2 . i 

L ' 77T7 E-2 + /, . _n!- 1 1! + ,, . _Nl 2 ! + * * * 


(d+a) 


(d+a) 


(d+a)' 


_i_ r e m,l-3 , *11.1-2 _a_ , g M.i-i ai + 

> 3 L (d+a) 1 ^ (dda) L - 2 11 (d-ta) 1 - 1 21 


'M.l 


r P + -rl'+jJa,— . a_ - JLl£. a— 
L P M,0 + (d+a) 1 ! f A ^ \2 2 ! 


2 

£U + 


(d+a)' 
2 


3 


d 


1+1 £ P M,C J- + p M ,l 2 -. + p m,2 


2 r 


n 3 

•Tj-f + ....] +....] (19) 
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Further ( 19 ) may be rewritten as 
a> i I*- 1 

'• - 1 5 , wtj ’ T , ?- 1 1 fr 


+ ^ L "' 1 q 2 , ] ,r P M,L-2 d 3 ^ 2 + P M 8 h-l_ d 1 " 2 q 

, 1,4-1 2! t ...i.ju _ T ,_9 . . t._i tT 


(d+a) 


(d+a) 


L-2 


(d+a) 


L— 1 1- 


+ 


T Xi***2 2 P,_ P 

. da . . ^ r -n . M.l a M,2 q2 

(d+a) 1 2! ] C M >0 11 (d+a) 2 21 


'Mxk. 


+ ....] + d [ P M>0 fr + P M,1 -2-j- + P. ~3 

3 


2T + X M,2 ^=7- +....]+....]• (20) 


collecting equal powers of q, equation (20) may be rewritten as 


co 

Pe= K 5~" 


r=0 (d+aFd 1 i! C C Pm ' L - 1 ' a+a 


( - <L ') p ~ 1 _L p ( -iL- 

V j ■- > + r TVr T.—O'.J 1 +• « 


■M,L-2 v d+a 


+ p ( + p 1 , _2-_ r p ( — ) L + p 

X M,1 ^ d+a J + X M,0 J + d 1 ! l x M,L ^ d+a ; + X M,1-1 M+a ; + 


+ p ] + t -.fl_ ) 2 i_r p ( JL. ) l+1 4.p r 

*** M,0 J 1 cL* 2! L ^M,L+1 ^ d+a ; + %,1^ d+a ' 


+ P M,0 ] ....+ ( -§-) 3 4,[ P H I+J ^ ( jfe ) 1+j_1 +. 


+ p. 


M,L+j-2 v d+a 


( -3. ■jL+j~2 

V r) . „ / + • • • • 


+ P. 


M,0 ] +* . , » . J 


( 21 ) 


Taking together the last L terms from each of the braced terms, 
one obtains 
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Pe = K 


js: 


±=0 (d+aA 1 x! C C Sm ' 1 - 1 ' r+a 


( ^ J 1 " 1 * P 


+ ^ 8X P (d/d ) + ( ^ ) yt E M f L ^ <d+a ^ 


•V 


d ' 1 


M»L-2 V d+a 
d \L 


( T 7„) L ~ 2 + 


+ ( “d" jZ h t E m, 1+1 ( dfe ) I+1 + P M,L ( d+£ )L ^ + ' 


+ < -i-) J b [ E 


M,L +d -l ( 5^ )W + P K(1+ ._ 2 ( ^)^- 2 


+....+ P K L ( ) L ]+....] 


(22) 


a 


Since 


? T, n - v _ ( L+k-1 ) [ <|> (L+k,Lj x) 

i+i,X ± , - (x^ITIk! 


(23) 


Where c£> (L,i,x) is 
series expansion 

<Js(L,i;x) 


Confluent Hypergeometric function with 



(i.) 


TTJ 


m 


n 



(24) 


in which (P) m = P (P+l) • . . . (p+m-l)for any realp 
Letting x = (g-|~ ) and summing each term of (22) over index 
i and using (23), one obtains 


Pe = Zsxziatel . . . [4(21-1, 1,;*) 

d'^fd+a) 1 ' 


d nL-1 ( 2L-2 ) * 

d+a ' (L-l)!(L-l) ! + 


4>(2L-2,Ljx) 



L-2 



+ 


+ <f>(L,L;x) ] 
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+ 


K 


[ [ esp( q/d )-l] ( ^_) I 4(2I„L;x) 


(d+a) L d L 

+ [ exp (q/d)“ -^-9 — ( — 1 ] ( ) <J>( 2H+1 , Ljx) 3 

(25) 

Substituting for K from equation (ll)one obtains from equation(25) 


L-l 


Pe= (jt-)! exp(-p/a) [ I4(I + i fIiX ) ( j^) 1 {i^> 


i=0 


1=0 


^>(21+1, ijx) ( {UxyrfpiaJ it i- fia ^ a ^-exp(-q/a)H 

( 26 ) 


Equation (26) may be used to compute the probability of error. 

As may be seen from equations (37) and (46) the ratio of (m+l)st 
and nth team of the series is less';than 


( 1 + 



) ( -2— ) 
; 1 d+a ' 



21-tm \ 

L-tia ' 


Hence the rate of convergence of the series in 
equation (26) is at least same as that of exponential series 
with argument (1+ — ^ Hence in the situation where 
(1+ -~-) ( -“L- )< 1 which is the case in most of the systems 
operating over Rician channels of practical interest, only 
a few terms will suffice. Each term can be written as a finite 
summation as given by equation (Al) of the appendix A, namely, 


(l-fm,L;x)= exp(x) 


m 

& 


x 1 ( L-l ) ! j! 

( L+m-1 ) I i !(m~i ) ! 


L>0 
m > 0 


(Al) 
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As an illustration, equation (26) has been used to 
calcinate Pe for coherent PSK modulation. 

Since (1+ is high for high ratio of 

(y g /Yd) so as to require large number of terms, a reasonably 
high ratio has been taken to calculate from (26) 

4_ 

Let L=4 • Y d = 1 , ^s-i = 16 

Prom (80) and (81) one obtains 

q/d = 2.34, x=2, (^)= .85 

To calculate Pe from (26) only seven terms of infinite 
series are required to yeild exact probability of error 
Pe= 3.74x 10“ 6 

If however L=4, Y^=4,^Y c; . 1 -=24 
then q/d = .32, x=. 3 

only two terms of infinite series are required to 
yeild Pe given by 

Pe * 1.77 x 10~ 6 


Equation (26) may be written in some what different forms by 
using the equality 


exp(q/d)- r 

J=0 31 


- 8 - > 
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to obtain 

P e ~ ^ d4a ) le2: P(“p/ a ) [ <f(L+i,ljx) (j^) ;L “[L-i)li! 

laO 

CO P T T . 

+ (-£» exp (-q/d) y ^ ) 1 i|-(l,i+2,q/d) • X 

4> (2L+i,L,x) ] (27) 


By Eu_nur bs transformation [32] 

<|>(l,i+2, q/d) = exp(q/d) <|>(i+l,i+2, -q/d) (28) 

Hence 

Pe " ^d-fa ^ exp(-p/a) [ jf^<f>(L+i ,L,x) (^ a ) =L 
00 P • 

+ ( -f-) £ TliTF ( d+a — ) i 4 5 (i-+l»i+2 7 -q/d) 4(21+1, L,x)] 

1-0 ' ' (29) 


2 . 3 Large SHE Assymptote s 

Substituting (24) in equation (26) and rearranging 
terms or directly from equation (21), one obtains 


Pe = ( d^^ L ex P(“P/ a ) ex P (~l/ d ) 


LsM: 


3=0 


3 1 


oo i 

X fr 


M=L+i 


L+.i-l 

X 

&0 


d xk 


x X p M t k ^ d + a ^ 


(30) 


Writing (~^ ) as [l- ] and using binomial expansion. 


a 


one obtains 
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t « oo 

Pe= ^ ezp e^p(~q/d) .( g/d ) J 


0=0 


0 ! 


1=0 
M=L+i 


x 
1 ! 


L+j— 1 

Z 

k=0 


( i ( M+L+j-1 ) ( M+k-1 ) ! / a \k ( , zi \ 
K ~ XJ ( M +k ) (M-l) !k l^d+a ' K:>±} 


How 


x__ (M+L+3-1) ( M+k-1) 1 ( 2L+.j-l) i(L+k-l) ! Y 

fz 0 i ! ( M + k ) (M-l) Ik !” (L+k) l(L-l) !k !(l+;j-k-l) ! 


(L+k,21+;j*L,L+k+l;x) (32) 

where 2 P 2 ^ #Y >6 ;x) is generalized Hypergeometric 
function [32] and 


«> W ).( P ), 
2 P 2( <?M3,Y,<5,x) - 2^ t7^T3TT 


x 


i ! 


(33) 


Changing the order of integration in equation (31) and using , f) 
(32), one obtains » 

Pe= «p(-p/a, exp(-q/d) 


X- TS^ k -i) ' ( - 1)k c afe 5 V 2 (i+k,2i +3i i,i + k + i } x) 

(34) 


It may be seen that for the slow fading situation 
( £~) represents a quantity which decreases as the signal to 
noise ratio increases. 
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Hence for large signal to noise ratio conditions, 

)«1, retaining only the lowest power of ( ) in 

equation (34), one obtains 

Pe= ex P(“P/ a ) exp(-q/d) ^ . (q/d) 3 ' (2L+j-l )' w 

a+a j=0 j! ? L+j-l)! 1 I A 

$(2L+j, L+l,x) (35) 

2.4 Upper and lower Bounds on Pe s 

An upper and a lower bound is found for the probability 
of error which is valid for the entire range of signal to noise 
ratio. Let the first tern of (26) be denoted by Pe^ and the 
second term by Peg, then 

Pe = Pe-^ + Peg 

Now Pe^ is upper and lower bounded by a single <|> 

function. As shown in the appendix, we have 

<f>(l+i,l,x) = (l+i— 1, l,x) + cj5 ( L+i , 1+1, x ) i > 1 (35) 

L > 1 

Further as may be proved by term by term subtraction 
of the series expansions, one obtains 

Cp (l+i-l,L,x)>d>(L+i,L+l,x) (36) 


Using the above inequality one gets from equation (35) 
cj) (L+i,L,x) < (1+ g“)^(P+i“l, l,x) (37) 
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Since (p(L,L,x) = 21 _Jl! = exD ( x ) 

X m=0 TlX m! exp w 

m 

by repeated application of equation (37) one obtains 
<|>(L+i,L,x) <C (1+ ) i <|) (1,1, x) 

or <|)(l+i,l,x) £l+ ) i exp (x) (38) 

Similarly substituting equation (36) in (35) again, one obtains 
<f> ( L+i , L , x ) > (1+ |j--)<J> (l+i,l+l,x) (39) 

By repeated application of the above inequality one gets 

<f>(l+i,l,x)> (1+ |-) (1+ x^i )••••( 1+ xfzi”) $ ( L+i ^ L+i » x ) (40) 


or we have 


and 


<|> (l+i,l,x) ^ <p(l+i,l+i,x) = exp (x) 
$>(l+i,l,x) > exp (x) 


(41) 


Substitution of (38) and (40) in equation (26) yeild the upper 
and lower bounds respectively which are given by 


-p t- sfe ] c K 


1-1 


k=0 


d+a 


4 ) k (i+^-) k + 


< 1+ rr )L < dfs ' L 


co p / d 


) Jj exp (~q/d) 21 x l»l+i v d+a 

i=0 




where 


e ± = exp ( q/d ) 


i 

S> 


laZU 1 

i t 


i ' L' 

(42) 

(43) 
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, , 1-1 
Pe ^ ^d 5e"^ exp afd+a ) ^ ^ L,k 

^d4a“) L ex P(“l/ d ) ^ P l,l+i ^d+a"^ e i 

i=0 

, i-i 

» where h(l,i,x) = JJ (1+ — ) 


( afe") 1 * 1 + 

h(l,l4i,x) ] (44) 

k * 

(45) 


Further it may easily be shown that 



(46) 


Substitution of (37) and (46) in the second term of equation 
(26) yeilds 


£e 2 4 ( 4+a“ ^ ex I’(~s/ a ) 


Jr (1 + -| )i <$>(2L,L,s) (^) 1+1 x 
f 2L+X-X ) ! a 1 

eimiST ! d i (1+1) , 


or 

Pe 2 ( a4a ~ ) L exp (-p/a) <J>(21,l,x) ( — ) L 


CO 

Z d + 
1=0 



or 

Pe 2 ^ 


X ( 


a )i 
d+a 

s 


1 72140 - 1 ) ■ 

T1 (1-1 J !( 1+i ) ! 


d^) 1 ^d4¥”^ L ezp ^“ p / a ^ 4>(21,1 ,x) cf> (21,14-1, (~^) (l+~f] 



x/ 21 ! 

A 1 !( 1-1 ) l 

(47) 

Further we have 

e > ( -a_) m+1 i . 

e m ^ ^ d ' (m+1) ! 

(48) 
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Substitution of the first inequality of (41) and (48) in second 
term of equation (26) yeilds 

Pe 2 > ( d+a“^ L ( dfe^ exp exp C— q./d> cf>(21,l,x) . X 



d \m 
d+a ' 


( -£-) m+1 


1 ( 21+m-l ) l 

Tm+1) I (1-1 ) i(L-tm) ! 


or 

Pe 2^ ( dfe )L ( dfe“ )L exp ( “ p/a) exp(-q/d) $(21,I f x) 


\ 


X [$(21-1,1, gqg-) -1 ] (49) 


Either of the inequalities (42) or (47) may be used 
to obtain the upper bound. Similiarly one of (44) and (49) will 
provide the lower bound. Depending upon the values of x and 
q/d one or the other may be tighter. When these inequalities 
are used for coherent PSE modulation over slow fading pic ian 
channel for y d and y g used earlier to calculate Pe from (26), 
one obtains 

(i) 1=4 , y d =i , ^ y si = 16 

from inequality (42) Pe^ 5.1x10"^ 

from inequality (47) Pe^-=2.72x 10"“^ 

Similiarly from inequalities (44) and (49) one 
obtains the lower bounds given by 
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?e L = 2.36 x 1Q“ 6 

£ 

Pe^ = 1,0 x 10“ respectively 

—6 

whereas exact Pe=3. 74x10“ 

4 

(ii) 1=4 » 

a i=l 

from inequality (42) 
from inequality (47) 

Similiarly from inequalities (44) and (49) one 
obtains the lower bounds given by 

Pe L = 1.76 x 10“ 6 

and Pe^ = 1.74x 10“ respectively 

whereas exact Pe = 1.77x 10 ^ 

Hence the bounds given by (42), (44) and (49) are 
quite tight especially for low values of ( Y g /Yg)« 

Figure (3-10) plot the exact probability of error and various 
upper and lower bounds on Pe for coherent PSK signalling for 
various orders of diversity. 


Pe n = 1.78xl0“ 6 
Pe n = 3.64 x 10“ 6 
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CHAPTER 3 

EQUAL DIFFUSED COMPONENTS 
FSK, DPSK AND PSK SIGNALLING 

3.1 In this chapter the probability of error expressions 
are specialized to various slowly fading situations. In parti- 
cular the error expressions are obtained for FSK, DPSK and 
coherent PSK signalling for slowly fading Rician channel with 
diversity reception. The results are further specialized to the 
situations where the specular component is zero and the received 
signal has Rayleigh distributed envelope. These results are 
shown to conforn to the ones derived earlier in the literature. 

3.2 FSK Signalling ; 

The variables a^, d^P^ arL( ^ 9^ c °ulcl be evaluated by 
appropriate substitutions in equations (6) to (9). However in 
this subsection these are derived in a direct manner assuming 
from the beginning that the channel is slowly fading. 

In FSK signalling scheme the received signal over 
the ith diversity branch is given by 

v i (t)=V i Sin w ± t+ R i (t)Sin[w_.t+ cp^t ) ]+n ± (t ) , (50) 

O^t^T, i=l,2,...,L , j=l,2 

2 % 

where w^ and w 2 are frequencies seperated by — 

■\A sin w^t represents the specular component 

R^t) Sin [w -t + <p. (t) ] is the diffused component of 
3 received signal. 
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Each of n^(t) is an additive narroxtf hand Gaussian process 

p 

assumed to he independent with mean power a . 

For slowly fading channels v^t) may he rewritten as 


v ± (t) = R ± Sin [w.t+ Q.] 


p(Ki_)= 


R- has Rician probability density function given hy 

2R. RT + V. 2R. 

exp [ - — ] I n ( 0 (51) 


p? V 2 

* X X 


el T ? 


= 0 


where 


( 1 + 




p ! 2 v? 

V 1 1 

o' 2 

— r ~ T 

Y ± /a 


0 p? V 
' x i 


) 


Ih< o 


p 

(p. 2 being the power ratio ~ 1 ^’^ le ( 3-i^f use( 3- ^.to 

v| 

specular component of the received signal. 


Iq ( ) is the modified Bessel function of order zero. 

The multichannel receiver for FSK signalling would 
be considered to consist of a matched filter and square law 
envelope detector for each of the two frequencies for each 
diversity channel and combiner of the detector outputs 
corresponding to various diversity channels; The decision rule 
is to choose the symbol whose corresponding combiner output 
is larger. 

The output of the combiner corresponding to the 
actually transmitted signal is given by 
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y = Z 


(52) 


2=1 


where y^ is the output of square-law detector 
corresponding to ith diversity branch. 

The density function of each of the random variable 
y^ is given by 

1 ' - ' »2 >, , y _i 

,2 

i ' i 


2 V> ± 


=„-Ttr- ex P (- 1 / ?i 5 ox p VsVPb ( aSZ — * 


p ‘ : Vf " ' 1 " pf vf ’ u P?V 

\ i i Nil Nil 

y i>° 


0 


ly (53) 

y i<0 


Assuming equal diffused components and uncorrelated 
fading over L diversity paths and defining 


a. = 


1 


1 P? V? 

' 1 X 


, b,= ( ) 2 f c i= —X- ezp(-l/ef ) 

Ni i 


P? Y. 

V X 1 


(54) 


(53) may be rewritten as 

P ( y n *)= C -; exp(-a ± y i ) I Q ) 


i" i - 


(55) 


For the case when all a^ are equal to a , from 
equations (52) and (53) the two sided Laplace transform of 
p(y) is given by 


afl P(y) ] = CT7°i] 


i 




i=l 


( s+a ) 


L ex P [“TiV)] 


( 56 ) 


The output of other combiner will have only noise 


component given by 


n 


L 

SI 

fa. 1 
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Wliere the probability density function of each component 
n^ is given by 

1 O 

p(n ± )= - — ex P (-3 ^/ct ) 
a 

Assuming that the n^’s are independent and identically 
distributed random variables one obtains 


[ p(-n) ] = ( [ 


*-S4-l/ <3 C 


]- 


The variable w is given by 
w= y-n 

Since y and n are independent, the two sided Laplace 
transform for the probability density function of w is given by 


*(s) ^X[p(w)]= ( TT °i ) d L — i~ T exp [- 

i=l (s+a)^ 


Comparing equation (57) with (11) yeilds 

L 


Zb . 


s+a 


] - 
(-s+d)' 

(57) 


d = 


L 

Z \ = 

i=L 1 
1 


2 _ Y 
i=l 


si 


a 2 (l+y,) 2 


q=0 


a= 


'd • 
1 


(53) 


cy (1+Y d ) 


where y g ^ is the specular component of th^signal to noise 
power ratio for the ith channel given by — and y^ is 


the diffused component of signal to noise power ratio given 

. • p ! 2 v 2 

by ( >-V"' )• 
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When the signal to noise ratio of diffused components are also 
unequal, one obtains 


where 


p( s ) = 7T P,(s) 
k=l * 


F k (s) =; cC [ P ( ) ] = & 


exp[- 


Pi 


k 


( S 4 Q.Tr-) 


(s+ a k ) 


1 

(-s+d) 


where K k =d exp (-p k /a fc ) 

_ J_sk _ n 

Ptr - p / ,9 » 9-k ~ ® 

« < 1+ W 


(59) 


a k~ JL 


1 


d+W 


, d, = l/o‘ 


where y., is the signal to noise power ratio of the diffused 
die -> 2, r 2 


component over the kth diversity channel given by £ 


fk 


V 


) 


Substituting (58) in equation (26) yeilds 


5e = tdfe“ )L e*P(-P/*) gj 


or 

Pe= 


(2+r d y 


exp 


r ^-^si -j ( L+i-l ) r \i 

1+Y d J ^~ 0 (L-l) I i! ^ 24-y d ' ' 

Y 

X $ (I+i ’ 1 ’ Ti^TTs+yTT 5 


( 60 ) 


d' vt - ,T d- 

Pigure (11-12) plot the probability of error as calculated from 
equation (60) for various diversity orders. 
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3.3 DPSE Signalling : 

Por DPSK signalling the variables f and g appearing 
in equation (l) are zero and e equals 1. 
hence 


w k= C ?k + 2 k £ 


and matrix & is given by 
& = 

Por the case of slowly fading channels the complex 


0 1 
1 0 


random variables and Y^ are given by 




( u 1 + n x ) +i (u 2 +n 5 ) 


~k = ^ u l + +D ’ ( u 2 +n 4^ 


where u^’ s and n^'s are statistically independent 


normally distributed real random variables. 

P’ 2 V 2 

u 1 = K (Y w ) 


u. 


k' 2 

f Y 2 
'k k 


2 = W (0, s ) 


( 61 ) 


n ± = 0 (0, 


o / 


i= 1,2, 3, 4 


Hence matrix and vector S are given by 
2 


H 


k 


(V V. + a- 


>'2 
l k 

)' 2 y «^ 

k k 


r 2 

'k 

r 2 


> O O 

fir v k 


i 


'* 2 

i'k 






S = 


Y, 


V, 
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The characteristic roots of matrix H k G are 


obtained as below. 


H k G = 


e’ 2 v 2 

v k k 


P'2 Tr 2 2 

V k \ + a 


?k 2 V k + ^ V 2 


k k 


Now the characteristic equation of matrix H k G is given by 


H k G- I>J =0 


or 


or 


1 




k *k 
-'2 


K -A 


£ V 2 + o 2 


k k 


€ t 2 ■ - 2 1 


k k k 


+cr 


C 4 


= 0 


-2 A f^ 2 V 2 + (<?l 2 V 2 ) 2 - (^ k 2 V 2 +cr 2 ) 2 =0 


‘k k vv k k 
simplifying the above yeilds 




2 A P^ 2 V 2 -(2 P^. 2 T 2 o 2 + a 4 ) =0 


k k 


The roots are given by 

v . - 2 ^k ± ir [ 4 <i 2 ^) 2 + 4 < 2 ] 


e ; 2 ± ( kx + a 2 ) 


or 


A; 

A- 


-a 


K z v 2 « 2 
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= ( A + <? 2 ) (A- 2 ^ 2 vjA-cr 2 ) 


comparing the above with equation (6) yeilds 


1 

i o p p 

2 ?£ v J- kA 


(l+2y ,, ) 


a k = 


(62) 


Substituting a k' d k and matrices G- and in equation (7) 


one obtains 


— — 1,1 1 1 ( 63 ) 

(2 P k 2 ^|+^ 2 ) [l 1 


finally from equatipn (8) one obtains 


= 


(2 ?l 2 vl+a 2 ) 2 


•k k 


ff 2 (1+2 t^v) 2 


(64) 


When the signal to noise ratio in the diffused component 
over various diversity channels are equal, one has 


- A. Pi 


and 4 = = 0 


a 2 (1+2 y 4 ) 2 


( 65 ) 
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Further for this condition 
1 


a = a k = 


d k " 2 

a 


CT ( l+2y d ) 
1 


( 66 ) 


where a and d are the variables appearing in (ll) 
Substituting for a,d,p and q in equation (26) yeilds 
the probability of error for DPSZ signalling over the slow 
fading Rician channel. 


Pe 


PTr . 

r ■ ■ * mn 


1-1 l+2y 


d 


2 L (l+y d )- 




P l,k X : 


(I+k - 1 - (in, 1 ) (i + 2y d ) > (67) 


It nay be observed that (67) could be obtained from 
(60) by replacing y, by 2y, and y . by 2y . . Hence DPSK 

CL Cl SI SI 

signalling requires 3dB less signal power compared to FSK 
signalling to achieve the same probability of error. 


Though equation (67) itself is used for computation, 
it will be put in a different fora so as to derive an assymptote 
for the high SHR condition. Using (24) in (26) or directly from 
equation (21) one obtains by putting q=0, 

CO 

5e = < aV )L e =p(-s/ a > 1: 


i=0 
M=l+i 


i! 




d 


P ( 

L+i,j v d+a 




( 68 ) 
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Writing d - i _ a 

d +a ~ ~ d +a 

in (68) and using binomial expansion of each term and 
rearranging terms, the coefficient of (-l) m ) m in the 

second summation of equation (68) is given by 


(M-Ha-!) ! ( m) , (M-tm) * ( m+1) , 

' m !(n-l) ! ( m) (m+1) l(M-l ) ! ( m ) + 


+ 


(M+h-2) I 
( l-l)!( M-l ) 2 


(L-l) 
( m ) 


which may be rewritten as 


(M+m-1) ! , (M -tin) ! (m+1) 2 (M+L-2) ! (h-l) I 

m !(M-l) : (m+l) 2 (M-l) 2 ml! + + m !(l-l-m) » 


_ ( M-Hm-1 ) 2 r ( M+m-1 ) 2 (M -to) 2 , 

- m !(M-l) ! L (M+m-1) lO ! + TM+m^lJTV 


(M+h-2) 1 

(M+m-1) Y(li-l-m) !-* 


How using the identity [33] 



a+1) _ ( a-n) 
r+l) ~ ( r+1) 


(70) 


where n and r are any nonnegative integers and a is 
real. Substituting a, r and n by (M+L-2), (M+m-1) and 
( L-m-1 ) respectively, in identity (70). the coefficient of 
(-i) m ( is S iven 
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(M-m-1 ) ! ( M+L-l) f M+m-1) 

m !(M-l) ; ~ [(«-«“ > " {Mho ) ] 

Since the second term in the braces in ze r0 » the 
coefficient of (-1) 111 ( — ) m within the second summation 

in equation (68) is given by 


(M+m-1) 1 ( m+L-1 ) 

m t(M-l) T ( M+m ) 


Substitution of 'the above in equation (68) after 
writing binomial expansion of various terms, one obtains. 


Pe 


a 

d+a 



00 

exp (-p/a) JjT 
1=0 



( M+L— 1 ) 

Jpr 0 ( M+j ) 


X 


1-1 


M-iT 



(69) 


Now since 
2 E 2 


a> 

( a'»P;ir»6>x) = X 

i=0 


Wl <«± 
Ty7~T67— 



(70) 


One may write 

(PC 

y _xi ( M+L-l ) (M+.i-l) ! _ ( 2L-1 ) l _ v 

i~0 i ! ( M+j ) (M-l) !j T ( L+j ).{ 1-1 ) A 


2 P 2 (L+0,2LjI+j+l,L,x) (71)' 

Changing the order of summations in (69) and using 
equation (71) we obtain 
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Pe = 


( jfc-jW-p/a) g (-1)3( ^ 


1-1 




\s?) m 


X 2*2 (l'+d*2L;l+3+l, l,x) 


(72) 


3.4 High SER Case : 

Under large signal to noise ratio in diffused 
component y^, one has 

a 1 


and 


d+a 


a 


i+2y d 

1 


d+a l+4y 


d 


« 1 for PSE signalling 

« 1 for DPSE Signalling 


Hence retaining the term corresponding to lowest 

power of ( ) and neglecting higher power terms, ';the 

(72) yeilds 


3&J -) 1 exp (-p/a) tifetfi&ZL.™.,*) 


(73) 


Substituting the value of a,d and x from equations 
(58) and (66), one obtains for PSE Signalling 


3?e= ( exp 


[- ^ Ysi ] t (21,1+1, ^ 1 )( 2+y ) ) 

L 1+v. J T V T>— 1 ) i V-Lvy 


1+y, J 1 1(1-1)! 


a /v--r d . 

(74) 



and for DPSK signalling 


Pe = 


(2+2 y„y 


exp 


[_ 2 ^s± 


I ( 21-1 ) ! JLf OT T -j Zv 

si 

-> 11(1-1) ! T 2Lj L+1 5 ( l+2y J)( 1+Y 


l+2y 


d 


(75) 


5.5 Pure Rayleigh lading PSE and DPSK Signalling; 

When the specular component is absent on all the 

2 2 2 2 
diversity branches, we have o&, TP — 0 such that 

is finite. As may be seen from equations (58) and (65) for 

PSE and DPSE signals p=0. Substituting for p in equation (60) 

and (67) one obtains the probability of error expressions for 

slow fading Rayleigh channel. 

For PSE Signalling 

Pe = — I tjr 1 (iH-k-1) ! , X t Y d \k (76) 

(2+Y d ) L k=0 (1-1) !k! 1 2+y d ' 


This result is the same as obtained by Pierce [14] . 
Por DPSE Signalling 


Pe = T 

2 1 (l+T a ) L 


X Clrt-k-l) ! , •* ' +2T a 'k 

(1-1) !k 1 v 2+2y r " 


( 77 ) 


3.6 Coherent PSE Signalling : 

like DPSE signalling scheme in this also, the 


variable w^ is given by 
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hence G- = 


0 1 


1 0 


In the case of slowly fading channels the complex 
random -variables and are given by 


Uv = (u 1 +n 1 ) + j (u 9 +n Q ) 


2 2 ■ 


^ + ju 2 


where u^ and n^ are statistically independent 
normally distributed real random variables. 

„-.<v 

p ' 2 Y 2 

u 2 = B ( 0 , -~- 2 - ) 


^1 * n 2 — ^ ( 0 1 2 ) 

Hence the matrices and S are given by 


f ' 2 V 2 

k - k 


1 2 

(i ^ 

t 2 ^ 

k k 


Further 


k k k k 


e' 2 v 2 « 2 


F* Y i < 

jow the characteristics equation of the matrix H^G- is given by 

K° -ia| =o 
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or 


» ? 

0 Y 2 
v k v k 


A 


P 2 T 2 


?k V k _ A 


or 


or 


A 2 -2 A + ( C v£) 2 - 


k k 
’2 


k *k' 
2 


A 2 - 2 A (1 V 2 - P ‘ V 2 a 2 


0 


?k ^ (fe 2 > 


0 


= 0 


The roots of this equation are given by 

*, 2 _ 2 ill ffifl 2 ** C T I g2 ] 

JL <j) C- * • rmrm.nr t r. ». f,-,-ir-_- am. -~m, urn tmmam , 


or 


*i= <?’ 2 V 2 - ? t( fl% 2 ) 2 + ffv 2 * 2 ] 


k k 


‘k 'k 


A,„= ^ V 2 + V [( P„ 2 T 2 )2 + P t % 2 0 2 ] 


■k k 


k “k 


Comparing the above with equation (6) one obtains 


or 


a. 


a 2 [ V (Y aI+ Y * v ) + Y* v ] 


a, 


k 


= [ 


'dk T dk y 
^ (Y dk +Y dk } 


dk 
Y dk 


Y dk 


( 78 ) 


and d^. = 


1 r ^ Y dk + Y dk^ +Y dk 
•2 L 


Y dk 
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Substitution of a^, and matrices G in 
equation (7), yeilds 
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Finally from equation (8) one obtains 

Y sk 

Sk ' 4k ~ 

For the case when = Y^ for k=l,2, L 

one obtains 


(80) 


a k -a 


i k =a 


i r Y a , 

^2 L v . J 


1 r < Y a + Y a>- Y a n 

2 L v J 


k= 1,2, . 


y r 


and also q= p= p 


k=l 


k 


H Y S k 


2a 2 Y d V( Y 2 + yJ 


d 'd - 


(81) 


Substitution of a d,q and p from (81) into (ll) and 
(26) will yeild F(e) and the probability of error for the 
coherent PSK signalling over slowly fading Rician channel. 

The probability of error expression for this case 
may further be simplified by substituting q=p. 

Substitution of (24) in (26) and arranging the terms 
in powers of (l/d) or directly from equation (19) by substituting 
q=p one obtains 

CO CO 


pe= -s-t z -i z 

(d+a) r =1 d 


1 


r i^O (d+a) 1 i! 


°L_ ^M, L-r+k 

*sz 


*k 


k=0 (d+a) 


l-r + -k k ! 


M=l+i 


(82) 



46 


Collecting terms having the same pox^er of p y (82) may be 
rewritten as 


Pe= 


K 

(d+a) 1 





_h_ 

(d+a) 


L-r+3 


21+3-r-l)! 
L-lML+jirj! 0 U l 


(2L+j-r-l)! (2L+.i-r-l) I -i 

L '.(L+j-r-l) ! 1 !( j — 1 ) > + ' (l+j-r)!(L-l) ! 3! 0 ! J 


(83) 


Pe= — — — T 
(d+a) L 




(d+a) 


L-r+3 




1+3 -r) (1+3-1) 

k ) ( j-k ) 


(84) 


Using the identity 


Id M !*! ' ! 


(85) 


k=0 


for any integers 3,N^ and such that 

3 £ x 19 s 2 


One obtains 
E 


Pe= 


(d+a)‘ 


Pt-. 


fe d r (d+a) L " r j =0 (d+a)J 31 L ^* L+5 “ r 


( 86 ) 


By changing the order of summations and making use 
of the fact that Pjj m =0 for m < 0, equation (86) yeilds 
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Pe 


K 


d(d+a) 2L ~ 1 


CO . 

r p j 

3=0 (d+a) 3 


i l 


L+j 

X 

r=l 


' d +a \ r-1 
' d ' 


'L+a,i+;j-r 


(87) 


Substituting back the value of E from (6), the 
probability of error for coherent PSK Signalling over slowly- 
fading Rician channel is given by 

Ee= (d+a - ) L-1 (■a+a 3I, ex P (-P/a) exp (-p/d) f ' ~J~ ' X 

3=0 (d+a) J j ! 

ttP / d+a %r~l p (88) 

1_ K & } *1+3, l+j -r 

r=l 

This result is considerably simpler than the 
corresponding result obtained by Lindsey [26] which requires 
Q functions and various i'jp functions which have the arguments 
different from the ones appearing in this thesis and can not be 
expressed by a finite summation similiar to the one given 
in (Al). 
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CHAPTER 4 

UNEQUAL DIFFUSED COMPONENTS 

4<.l In the previous chapters it was assumed that the 

diffused signal components received over various diversity- 
channels are equal, though the specular component might he 
different. However, in many physical situations such as in 
angle diversity tropo scatter systems, the diffused components 
of signals received over various multipath channels are 
unequal. 

In the following we consider a situation where the 

1 diversity channels could be divided into N groups such that 

the channels belonging to the same group have equal diffused 

components but in general unequal specular components. 

If m- is the number of channels belonging to ith 
1 N 

group then we have m. = L. 

1=1 x 

The specific case where each group has only one 
channel i.e. the diffused components in all the 1 diversity 
channels are unequal is considered subsequently in section 4.3. 

JLi 

Now E(s) = ~fT F, (s) 

k=1 r Pk r !v 

exp ^s+a — 3 ex:p [ a -s 3 

and ? k (s) = E k k '(-s+djj.T^ 1 

as shown in equations (5) and (9), one obtains after little 
simplification that 
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r n r ^1 % 1 
exp [ Ti^7) + -- + (^) 3 exp[ Ta^) + -” + (d^ s ) ] 


P(s) = E 


( s+a^) 21 !. . . , ( s+ajj) m F 


1 




(d 1 -s) m l...(d N -s) m ir 


(89) 


where E = J~[’ 


E. 


i=l 


m. 

>i= ^ P- - 

t=l ^ 


, q ± = 


m i 

iFi 


L i ,3 


(90) 


where p. - and q. . refer to the variables p and q 

X , J X f J 

for the jth channel of ith group. 


4o2 Derivation of Probability of error for 
P SIC/ DISK Signalling ! ' 

The probability of error Pe is given by equation(4) 


if 


e. 



p~(w) dw 


(4) 


and p“(w) is calculated by finding the sum of residues of F(s) 
in the right half complex plane like in chapter 2. 

In the following, first the expression for the 
probability of error is derived for the case of slow fading 
Rician channel with FSE or DPSK signalling. Under this condition 
as shown in equations (59) and (64), we have 


4 X =0, P ± 


4l =o, Pl 


sx 


1+Y di 

2Y si 


l+2y 


di 


for FSE 


for DPSE 
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Hence equation (90) yeilds 


q\ =0 for both FSK and PSK 


Y 


si 


and p i = ""57 77 — 75 

° a+Y dl ) 


(91) 


^si 


p± a 2 (l+2Y dl ) 2 


m . 


where = V 


si 


for FSK 


for DPSK 

i=l,2,...,N 


Y . and denotes the total 
' si; 


3=1 3 ^ 

specular component power in all the channels of ith group to 

noise power ratio in any of the channels. 

Further d.=d = -i- for i=l,2,...,N 
1 a 


Denoting Fe(s) by __ 

p-j^ Pjj ^ 

Pe(s) = exp [ T5 ^ i ) +..••+ Xi«„) ] 


(92) 


In this case, therefore, the equation (89) may be rewritten as 


]? ( s )- K /'/i I ' e ^ / . \m, 

(d-s) 


(93) 


(s+a 1 ) m l....(s+a iT ) m H 

Since Fe(s) is analytic at s=d, a Taylor series 
expansion can be written for Fe(s) around s=d as 

co 

1_ ___ 

i=0 ±! 


(94) 

Fe(s)= JP L S z^ LL- Fe^ i ^(s=d) for all s such that 


ld-sj < a min 


where a min =®ain(a 1 , . ,a M ) 


'IF 
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where 3?i i ^(s=d) = 

ds " Ld 

letting f(s) £ [ T h^y + ....+ T |^- ) ] 

we haTe -gsia)= f (Dp e , = [ [ f (l)]2 +f (2)]P eCs ) 

ds^ 


d^FeCs) 

ds^ 


where 


and 


-4-^ = C [ f^] 3 + 3f f(3)]y e ( g ) 


= [ [f (l) ] 4 +6[f (1) ] 2 f (2) +4 f (1) f (3) + 3[f (2) ] 2 +f 


, 1 = 1 , 2 ,.... 

IT — 

= (-D k W X P 1 

s=d 1=1 (d+a i ) k+1 


.(i) 


jA 


ds‘ 


■ (k) 


(d) 




d k f 


ds 


k 


(95) 

(4)j 


(96) 


The rational function in s can be factored by partial fraction 
"tpansion to yeild 


3? (s) = 
z 


1 

( s+a 1 ) m l . . . . ( s+ajj-) m N 


where 


N 


m . 

r 1 


In 


( s+a ± ) 1 


A i,3 ~ 


i=l j=l 

THPP zspr } 


ds i 


\ ' 


(97) 

*98) 



' s CENi.t/\L LSBRAR 
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letting 


) i 

(d~s) L (s+a ± ) 3 

Fe( s) 

(99) 

one obtains 

?(s) = 

N EL. 

X X i ( 

i=l j=l 1,3 

s) 

(100) 

Now letting 

p i,3 (w 

>n of" 1 [ ] 



We have 

p'(w) = 

N m. 

X X Pi -j M 
1=1 o=i ,3 



where 

pr i(w) 

± 9 J 

= 0 w ^ 0 





= p(w) w^O 



One obtains 

from equation (4), the probability of 

error 

Pe = 

N m. 

x x 1 

i=l 0=1 

0& 

f Pi j (-w)dw 

J 0 


(101) 

Now 





p7 ,-(w)= the 
- 1 - 1 j 

coefficient of in the 

Laurent 

series 

expansion of 

F. . 
if 3 

sw 

e 


(102) 


using the identity 


1 

( s+a^ ) 3 


= (d+a ± )~^[l~ 


l.-s -i-j 
d+a^ - 1 


(105) 


Substitution of (94) and (105) in equation (99) yeilds 
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F, , e sw = A, „ 

1 9 J 


(d+a. )“' 3 co , d-s -] . , 

1~ E Fe(d) [ ^j,r ^d+a. ' -J • X 


XD (d-s) 


r=0 


t k ** c L < d -> r %f3 eWd ^ 


k=0 


r=0 


where 


) k 

Fe(d) 


r._ 4 (-!)”• d^Fe(d ) 


ds' 


s=d 


(105) 


Multiplying out the first two series in equation 
(104) and rearranging terms one obtains 


. , A. , E Fe(d) cc co . / k 

, et sw = -j- [ I f n r (a-s) r ] [XI (a-s) k ifh]e wa 

" L * J /* r=0 1 9 J 9 ™ k=0 K * 


(d+a^) ^ (d-s) 


k 

(106) 


where P i>;j(0 - % =1, P lfj(1 = [ flSsp *j,i + TT ^ 


Oo P-. 


ax. 


ft - r -i .2 . r.i.l . *2 1 

X,2 " [ (d+a,) 2 1! -2 ! ] 


and in general 


y Jj xikzg. 

r^O (d+a^) r (k-r) S 


i — 1 f 2 ^ * • • « y E 

J 1 m_-^ 

k=l 1—1 ) 


(107) 


How, from equations (102) and (106) one obtains 

-A. , E Fe(d) 1-1 j i v 

p- ( w )= — exp (wd) J2 p . k (-w) ~ ~ 

1,3 (d+a, ) J k=0 1 » 3>1E (l-l-k) » 


(108) 
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CD 1 

Using the result f w exp(-wd) 

J 0 k ! aw ~ 


1 


. k+1 


one obtains from ( 101 ) b 7 substituting (108), 


E 

Pe = 


“i 

t. 


A. .K Pe(d) 

J — . ... T 


i=l j=l (d+a ± ) J 


1-1 

r 

k =0 




d 


L-k 


(109) 


For FSK signalling by substituting (91) into (92) one obtains 
for s=d. 

If K 

Sl 1 (no) 

di J 


^ f i Ti^7TU^-7 J 


'di' 


Further from equation(58) and (91) one obtains 


K= 


1 

n t 

±=1 x 


E 


E 


( -^) 2 I exp[- Z (fe— ) ] n — n_ 

a 1^1 ±+Y di • i=l(l+y di ) m i 


( 111 ) 


Substitution of (91) in equation (96) ye-ilds 
f (j) (d)= ( cr 2 ) ^ ( - 1 ) ^ j! ]T T sk ( 1 +Y d k ):j 1 


k=l °' a ' r ^ , i j +1 


Letting 


T 4 


■( j) 


( 2 +Y dk } 


(d) 


(o 2 ) J (-D j 


( 112 ) 


one obtains 


E 


l - 3! t y 


k=l 


Sk 


^ 1 +Y dk^ ^ 


( 2+Y dk> 3 ' +1 


(115) 
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Normalizing 8. . , and R. by defining new variables 

1 ^ J > ix 1 


A. P -i m* 


= and R. = — 

C (a 2 ) k 1 (<? 2 ) ; 


'i,j,k , 2sk — -I” ,2,1 


(114) 


One obtains from equations (95) , (96),(105), (107) , (113)and (114), 


„ , = 1 p 5^_ ( JZ&L f 

3 ,r (k~r) I 1 2+y, . ; 

r=0 ai 


(115) 


^0 1? R 1 S'!* R 2~ ^ ?1 + ^2^ R 3“^ fl + 5 ^i ?2 + ?^ 


R 4~ ^1 + 6 .fl 5 2 + 4 fl §J + 2 ?2 + ?4 ^ 


( 116 ) 


The R.'s required in computation of Pe with upto 

J 

quadruple diversity are listed in (116). Others could of course 
be derived from equations (114) , (113) ana (105). 

The coefficients A. . are calculated as given belows 

1 f J 


1,m i ( s+a^ ) m l . . . ( s+a^^ ) m i-l ( s+a i+1 ^ m i+i; : ; ( s+ajj) n IT 


s=-a. 

i 


rU 


which after substituting from equation (58) becomes 

A. - (a 2 ) 1 -"! fl [ fr-tW fr 

1 ' m i r=l ''a- 'in 


(ii7) 
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Further 


d 

' ds 

[ ( s+a ± ) m i 

I 1 

CO 

• » m i 

IT 

X . 

— n 

r 

V- 

r=l 

(a -a. ) 


r x' 




or 


A. 

i ,m^-l 


A i ,m. 


IT 


Letting 


3. .= 


n 

1 a 2 X 

r~l 


A iil. 


± ,j (0 2 } L^ 


- n r (l ^di )(l ^dr ) 
^ Y di" Y dr^ 


(118) 


for eu y- 1 
i=l,2,..,lT 


(119) 


One obtains 


B i,m, 


* r U +Y^Hi+Yfa.) m 

I I ( v Z v ) •* ^ 

^ lY di Y dr' 

r^i 


B. = B l. m i 

i , 321 . — 1 1 • 


IT 


-mJl+Y^)(l+Y dr ) 


r=l ^ Y di“ Y dr^ 


» m i > 1 


and in general 

t3 1 d m i~ J ‘ r /■ , 1 

i,j Tii-J ) '■ ds m i-3 t( 


1+ Yai _):i - Z<S> 11 


(120) 


X+Y di 


where F (s) is obtained by replacing a. by (t t t ; ) for 

£i J J ^ ,| 

0=1, .. ..IT, in F z (s). 

Hence substitution of (110), (ill), (114) and (119) into 


equation (109) yeilds 
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H H - F m. 

Pe = f f~! — i— jjesp [- Z yjg— y] II 1 B ± . X .. 

’ ' 12y ai ) i=l 3=1 1 ” ! A 


1=1 (1+Y dl ) 1 


1=1 ri ^- 


(4^-) 3 ^ "i.J.k 

^ + ’di k=0 


( 121 ) 


where r\. . , are calculated from equations (I15)»(ll6) and . 

(113) and B_. . from equation (120) 

X 9 J 

Expression for the probability of error in case of DESK 
signalling is obtained by replacing Y gi by 2 Y si and Y^p "by 2 ^di 

in equation (121) as in the case of equations (67)and (60). 

Figure (13-16)plot the probability of error calculated from (121) 
for dual space dual angle diversity troposcatter system where the 
signal strength over various paths are unequal. The plots are 
given a range of-Y^. 

4 . 3 Ho Two channels having Equal Diffused Signal Components; 

The general expression (121) could be simplified 
considerably in the case where no two are equal. This 
situation for example corresponds to the case of triple angle 
diversity with no space diversity on troposcatter link, in 
which the beam at grazing horizon receives maximum power and 
the one with maximum elevation angle has lowest received 
signal power [30], 
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In this case we have, 

a i ^ a k for 1 ^ k J i , k=l , 2 , , . . , 1 
N = L 

n^= 1 for i=l,2,...Ii 

B. • = 0 for j > 1 
± t J 



In this case the expression for Pe after substituting (114)'}! 
(119), (122) and (12J) into (121) becomes, 



(124) 
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After a little simplification of (124), one obtains 


Pe = exp [ - ST 


Y si L 


■ 4S ~— . 
i=l Q1 




l+2y 


1 ) C TT ] X 


i=i di k=i di dk 

kj4i 


L-l 


[ Z r H 


r=0 


i,l,r 


(125) 


Wh 


ere the ri . , are obtained from equation (115) as 
i,-L,r 


^1,1,0“ 1, ^1,1,1 ~ f R 0 ( 


1+Y di ^ . R i 


2+Y di 


•) + 


1! 


] 


and in general 


^i,!, 3 


1+T^-; R-: 

r -n ( > 4. ri , 

L n i,l,j-l ' 2+y di ; + 


] (126) 


3 , i — 1 , 2, ....,1 


where R' s are obtained from equations (116) and (127) 
where (127) given below is obtained from (115) by substituting 
R=1 and Y si by Y si » i.e. 



J! 


1 

21 

k=l 


sk 


t^Vak ) 3 -1 

t 2 + v ak ) i+1 


(127) 


Pigure (17-20) plot the probability of error as calculated from 
(125) for the triple angle troposcatter systems for the specific 
case of no specular component • 
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4.4 ' Equal Diffused Components i 

Since tile expression (121) is applicable to quite 
general situations, some of the results of chapter 3 could be 

derived as its special cases. 

Hence the result of equation, (60) follows from 
equation (121) when all the channels have equal diffused 
components as must be the case. 

In this case we have 


N=1 , m^= 1 


B. •= 1 for i=l, j=l 

- 1 - 9 J 

= 0 otherwise 


(128) 


Substitution of (128) in equation (121) yeilds 

1-1 


Pe = 

or 

Pe = 


d+Tai) 


exp [ - 


1 

r 


Y 


si 1+y 


Y 


dl 


] ( 


2+y 


ai_\i J" „ 


dl 


•)- 


k=0 


1,1, k 


si 


t. ex P C*"\j.2y -5 

(1+2 Y dl ) Y dl 


1-1 

^1 , 1 , k 

k=0 


(129) 


Prom equation (113) by substituting N=1 and letting y^j_ 


f^qual to Ypj » one obtains , 

... ( in *)*- 1 

3! Y s ( 0 ^, ^'+1 


where 


i 2+Yd) 

«*■ A y 1 

T s - X Y 
i=l 


(130) 


si 
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Further from equation (Hi), one obtains 


„ 5T (L+r-1) 1 R k~r 1+Y d , 

^l,L,k ' J= 0 (L-i')tr! Ik^FTT <■ 2+y fl ; 


(131) 


Substituting for L ^ from (131) into equation (129) 
one obtains 

L_ 

Y 


Pe 


(2+r d )‘ 


exp [• 


, si j; 1 y (L+r-l) ! ( 1+ Y d ) r R k „ r 

a —3 kto ho (2+Y , d) r irar 


l+2y, 


(132) 


Letting R. = 

v> 


R i 


5 2+Y d } 


and x _ (1 ® )(2 j 


(133) 


Substitution of (130) and (133) in equation (116) results in. 


Rq = 1 , = x 1 ! 


Rg = 2 lx + (1 ! ) 2 x 2 = x 2 ++2x 
R^ = 3lx + 3x 2 21+ ( m) 3 x 3 

3 2 

= x + 6x + 6x 

R 4 = 41 x + 3( 2!) 2 x 2 +41 13 lx 2 + 6 ( RJ) 2 2 x 3 + ( U) 4 x 4 
= x 4 + 12 x 3 + 36x 2 + 24x 


(134) 
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Substitution of (50) in equation (1J2) results in 
the probability of error expression, 


1-1 


Pe = 


(2+Y *)' 


exp [ - 


1+y. 


l+2y 


cl 


i / “'''a xk ( L+k-1 ) ! w 

J < 2^7 ~ } tmyri i X 


k=0 


■f* k- (L+r-1) • -p 

Z_ ( L+k-1 )! (k-r)! r > -k-r 
r=0 


(155) 


k 

Letting T,_ = k! (L+r-l) '• 

k XL+k-l71 Tk-r) 1 ! r l^k-r 


1 = 1 
o 


T 




1 L + Rq = 1+ -f = e“ x (L+l,L,x) 
2R ± Ro 


T 2 ^0 + Tl+i) + ~~i( L+l) 


2 

-i , 2x , x +2x 

~ 1 + TlS.T + “ l(l+i) 


= 1+ (L+1) + Li ( L+l ) “ e <i> (1 j+2,Ij,x ) 




5 a. 


^ ^2 S 3 

^0 + ( L+2 ) + Il+ 1) (L+2) + L ( L+l ) (L+2 ) 


= 1+ + rffen + r (i+i);u2T = e " x t < L+3 » L > x > 


and in general = e” x <§> (L+n,L,x) 


( 156 ) 
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Substituting equation (136) into (135) and substituting for 
x from (133) in the exponential term one obtains. 


Y 


Y< 


L-l 


Pe “(i^> [ " ] exp [ ‘(^T 


T" ( IrriC-l ) I 

^0 (i-l)!kT 


X <|)(I+k,I,-x) 


1+r a 'k 
2+Y a J 


or 


Pe 


(2+y a ) 


L-l 

T exp [- 5— — ] . 21 

L 1+Y d k=0 


L-fk-1 ) ■ / 
L-l) !k! 1 


1+Y d yk 
2+Y d ^ 


<^(L+k,L,x) 


(137) 


This expression is same as the one obtained in 
equation (60) 


4.5 ho specular Component and Unequal Diffused Components? 

One more specific case of the situation considered 
in equation (93) analyzed subsequently is slow fading Rayleigh 
channel with FSK or DPSK modulation (unequal signal strengths). 

For Rayleigh channel the specular component y . is 
zero for i=l,2,...L. Hence 

Y-, = 0 for i= 1,2,..., N 

oX 

Substitution of y . in equation (113) yeilds 

SI 

"§.= 0 for j=l,2,.... 

J 
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Hence from equation (116) we have 
Rq= 1, Rj£=0 for k / 0 


Substituting these values of R^ in equation (115) one obtains 


r , _ p ( ) k 

i,j s k 2+y di 


i=l, 2 , . * . ,U 

3 — 1 > 2 j # # # y m.. 

k=0 ,!,...,( 1-1 ) 


The coefficients B. . don't depend on y . and hence are given 

1) J sn_ 

by (120) as such. 

Substitution of m , in equation (124) results in 

1 » J * s - 


E 


E m. 


Pe 


= [ n 


] zi 1 B, , (■ 


1+ Y A "i 


i=l (1+Y di ) m i i=l 0=1 1,:5 2+Y di 


1-1 


k=0 


X ( 


1+Y di \k 


2+y 


di 


- ) J 


Letting 


B. . 

-it 0 


B. . 

l. -1 

B. m 
i m « 


for o ^ m. 


One obtains after substitution for B. and after a little 

jl f nij 


simplif ica t io n 


Pe= 


N m. 

Z_ B, 4 


1-1 


i=l 0=1 i,j (l +Y di) m i""^ k “° 


t— p /^‘ +Y di \k 

2_ ^0,k4+y r1 , ' 


di' 


'di 

( 138 ) 
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The probability of error for DPSK signalling is obtained by 
replacing y di by 2 ^ in (138) 


4.6 Rayleigh Channel (Quadratic Form Receiver) s 

In the following the function F(s) of (89) will be 
specialized to a situation where the specular components are 
zero over all the diversity branches. A general quadratic form 
receiver is considered from where the result for slow fading 
channel with coherent PSK modulation can be derived as a 
special case. The result for the specific case of PSK and 
UPSK modulation has already been derived in equation (138). 

For Rayleigh fading channel the mean of random 
variables and appearing in equation following (6) are zero. 

i.e. 


E [ 2fc] = 0 

k = E [ I k ] = 0 


hence 


^ k 


= 0 


Hence from equation (8) one obtains 


p ~l^_ — 0 f q^ — S f oh h — 1 ,2, . .. ,. 
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Substitution of and q^ in equation (89) results in 


PCs) = K 


(d 1 -s) m l (s+a 1 ) m l (djj-s) 21 !' (s+a K ) m N 


(139) 


The function F(s) has singularities at points s=-a^ , . . * , -a^ 

and s= d-^ , . . . , d^ 


Now as shown in equation (4) 

0 

Pe = | p“(w) dw 

-CQ 

Where from equation (10) 

p“(w) = - sum of residues of F(s) in the right half 
complex plane. 

N 

i.e. p“(w)=- residues of F(s) at 8=1^ (140) 

i=l 


P(s) being the rational function of s can be factored by 
partial franc t ion expansion to yeild 


N m. 

F(s) =K[I I AdL 
1-1 ‘ )-1 (s+a 1 ) ;] ' 


-f 


EL . 


r i 1 

i=l j=l 


D. 




(-sh<U ) 


^ ] (141) 
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Since each, of the term in the first summation is 
analytic at all the singularities d^, their contribution to 
the residue of f'(s) at the desired singularities are zero. 
Hence, one needs consider. . only the residue of the function 
F 2 (s) given by 

H m. t\ 

F ? (s) = K jr 5 1 -A t A (142) 

i=l j=l (-s+d^)' 3 

How the function ■ - — ■■ is analytic at all d, ,k= 1,2,...,H 

(-S4d i ) j K 

except when k=i. 


P^w) 


Hence from equation (140) one obtains 
H 3 D., , 

-E y residue of [ — — “ — ^ ] evaluated at s=d. 

i=i ( - S 4 d i )J 

(143) 


How 


the coefficient of ( g-- ) in the Laurent series expansion of 


— 7 e SW is given by 


(-s-Kl, ) 3 


e «L 

1 3 - 1)1 


N 


m * 

JL 


(-wd, ) j ~ 1 


Hence p~(w) = E X I 1 4 . wd. 

i=l j=l 1,3 (j-1) ! e 1 


Since 


0 

Pe= I p“(w) dw 

-CD 
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We have s using the result 

co 


{ _wL_ 
J 0 rS 


e^ iw - -jj- 


N 




Pe = K 1 T 

i=l 3=1 


D i.» j 

< d i>' 


(144) 


Now as one may observe from equation (140) , FgCs) 4s 
the Laplace transform of the function p~(w), i.e. by partial 
fraction exparsion the function F(s) has been split into two 
functions F^£s) and ^(s) wll ich are analytic in the right half 
and left half of the splane respectively. 

Hence from the mean value theorem of Laplace 
transform, we have 

co 

S P~( w ) dw = f 2 (s) | s=0 


From the above theorem, the probability of error Pe 
can be evaluated directly from (142) by substituting s equal 
to zero thus avoiding the evaluations of various residues and 
subsequent integration. 

The coefficients D. . are given by 

1 » 3 



a m i"3 

ds m i-3 


[ F(s)(-S4d i ) m i]| 


fe=d ± 


(145) 



Substituting for j =nu , (m^-1) and (m n . -2) one obtains 


F IT 

D± ’ m i = [ H T^r ] C T^T’^ 


m ± > 1 


r^i 


D. 


i,m^-l 


D i, a ± ■ t 52 -m k 




N 

X 


m 


k=l a i^k S. d r- a i 

rj^i 


h ^ ± >2 


D. 




^i,m. H— -m,, 

- 5 — — [ X [( ji-t 

k 


N IT 

-ra„ ni 


- t( )[ 1,-rl + z. ir=d-r] 

=1 i k r=l a r r=l v r i 7 

rj^i 


m k( m k +1 ) 


<W 


IT N F 

•] + [(3-4-) [ k Tertrr + X 7f-rr r 

k=l a i +d k r=l ^ a i + r ; r^l (a r i } 

k^i r^i 


m k (m k +1) 

(-d ± +d k ) 


■] ] 


which. after simplification yeilds 


D. 


D. 


i^i 

D i,m.-2= C< 


D i;m 


ilL.)2 + ^ m k( n k +1 ) _ ^ n k^ m k +1 ^ 


k=l 


k=i (a^ap 


kj4i 


( 146 ) 
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For the case of coherent PSK detection and slow 
fading condition, expressions for and d^ may he substituted 
from equation (78) i.e. 


E = 



i=l 


F 

= n w * 

i=l 

1 r 1 ^di* Y di^ ~ Y di 

2 L y -* 

cr Y di 

_ 1 j- ^^di^di^ + Y di j 

—2 Yaj 


(147) 


The coefficients a. ,d . and D. . could be normalized 

ii i , j 

2 

with respect to (l/a ) to yeild 


IT 


F m ± 


D. . 

Z3jJ. 


Pe = [ (a^ri ] 

i=l i=l j=l (d ± ) 3 


2 2 

where a^= a^c , d^ = d^cr 



and 3D. . is obtained fron 3D. .in (145) by replacing 

a^ and d^ by a^ and d^ respectively. 

3 - • 6 • 

[ P(s) (— s+d 1 ) n i ] i 

■ ’s=d^ 

F(s) is Obtained fron F ( s ) by replacing a^ by a^ 
and d^ by cL in (139) 

Figure (21 ) plots the probability of error for dual 
space dual angle diversity systems for the ca«e of 
coherent modulation with no specular component as calculated 
from equation following (147). Figures (22-25) plot the 
corresponding Pe for triple angle diversity troposcatter systems 
for various power ratios of the signals received over different 
antenna beams. 


D. 


■ . \H.-j .m. - 
, —1 ) i ** d a 


(13.-3): 


-1 m - 
ds 1 
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CHAPTER 5 
CONCLUSIONS 


The expressions for the probability of error for 
binary communication over selective fading Rician channels have 
been derived using the results from the theory of Laplace 
transform. It is shown that to calculate the probability of 
error Pe it suffices to calculate the probability density 
function of sufficient statistic w only for negative values of 
w. This is achieved by calculating the sum of the residues of 
the two sided Laplace transform P(s) of the p.d.f . of sufficient 
statistic w in the right half complex s~plane. The integral of 
this sum yeilds the required probability of error. In those cases 
where F(s) could be easily split into two f inactions E-^(s) and 
EgCs) which are analytic in the right half and left half plane- 
respect ively, even the evaluation of residues and subsequent . 
integration is dispensed with and Pe is obtained simply by 
evaluating E 2 ( s ) a f s=0. This will be the case for example when 
F(s) is a rational function as .shown in section 4.6. The 
expressions obtained for the various cases when the diffused 
components received over various diversity paths are equal, are 
considerably simpler than the ones obtained earlier which 
involve special functions like generalized Q functions requiring 
special computer routines. Though it has not been possible to 
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derive the one from the other, however, the probability of 
error calculated for a large number of cases from them is same. 

We have also considered a more general situation 
which in more complex in terms of analysis. In this situation 
the diffused components over various diversity paths are 
unequal, as may arise for example in angle diversity tropo scatter 
system. The probability of error expression is found for the 
case of Rician channel with FSK and DPSK modulation. Also 
evaluated is the probability of error for the general quadratic 
receiver operating over Rayleigh fading channel. The coherent 
PSK and FSK receivers being the special cases of the quadratic 
form receiver, the probability of error is evaluated for 
these hlso. 

This procedure will always result in the expression 
for probability of error containing only a finite number of 
terms whenever the function F(s) has only finite order 
singularities in either the left half or right half complex 
s-plane. In many situations which cannot be tackled analytically 
by method of convolution, this method could possibly be 
applied successfully. 



The probability of error expression derived for 
general quadratic receiver could be used to calculate Pe 'when 
there is inter symbol interference from adjacent symbols. 
Since in practice there is interf errence from only a few 
symbols the probability of error could be calculated for each 
possible sequence and then averaged out to yeild Pe for such 
a situation. 

As a ' byproduct a transformation of Confluent 
Hypergeometric function into a polynomical function has been 
derived when the arguments of Hypergeometric function have 
a certain form. 
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APPENDIX 

Here we give the proof of the following identity 
which is used in the evaluation of Pe. 

n 


4 (l+m,l,x) = exp (x) 5* x 

i=0 


i ( L—l ) • 
(L+i-1)! 


m i 

i! (n - i j ! 


(Al) 


for 1 ^ 1, m ^ 0 


Proof °° The proof is given by induction. 


oo 


Since ti)(i,j,x) = V 


X 


n 


n 


(A2) 


where (i) = i(i+l) . . . (i+m-l) 


By substituting i=j=L in (A2) r one obtains 


4>W.l.x) = Z -£r- = exp (x) 


n=0 


(A3) 


Further substituting i=L+l, j=L in (A2), yeilds 


4(1+1, i,x) = ]T 


Ob ( L-fi ) 


Try 


m=0 


m 


m x 


m . 


m ! 


or 


co 


4>(l+l,l,x) = 51 ( “Hr" ) 


111=0 


_x* 

m ! 


(A4) 
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Substracting (A3) from (A4) term by term, one gets 


<f>(L+l,L,x) -<J>(L,I,x) = [ 1+ “tQ 4 §r~ + * • • * ] 


JL. 

L e 


or 


<J>(I+l,L,x) = (1+ e x 


which proves identity (Al) for m=l and for L 1. 
Assume that (Al) is valid for m=n, i.e. 


4>(I,+n ( I, X ) = exp (x) x 1 5 X^ 3 ) , (A5) 

i=0 

L >1 


As can be verified by the series expansion and ten' 
by term subtraction we have 

<$>(L+i,L,x) - <f>(l+i-l,L,x) = — 4>(L+i,L+l,x) 

i ^1, L ? 1 
or 

<|>(L+i,L,x) = <|>(L+i~l,L,x)-f- — <|> (b+i,L+l,x) (A6) 

putting i=n+l in equation (A6) and using (A5)» 


one obtains 
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Hence (11) is also valid for m=(n+l). Since it 
has already been proved that (11) holds for n=0 and m=l 
it implies that the identity (11 ) holds for all non-negative 
integers, m. 
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